TAIL BEHAVIOUR OF THE AREA UNDER A RANDOM PROCESS, WITH APPLICATIONS 
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Abstract. The areas under workload process and under queuing process in a single server queue over 
the busy period have many applications not only in queuing theory but also in risk theory or percolation 
theory. We focus here on the tail behaviour of distribution of these two integrals. We present various open 
problems and conjectures, which are supported by partial results for some special cases. 
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1. Introduction 



In the past two decades there has been done enormous amount of work on asymptotics for queueing 
systems. Tail behaviour of steady state queueing process {Q{t),t > 0}, workload {W{t),t > 0} or 
busy period r in standard systems has been well-understood in both light and heavy tailed case. 
Surprisingly, however, very little is known on tail behaviour of integral functionals of the form 



where {X{t),t > 0} is a stochastic process (typically, X ~ Q or X = W), / is a deterministic func- 
tion and T is either r or deterministic (finite or infinite). Such integrals appear naturally in analysis 
of ATM. The reader is referred to references given in ||9l and ||22|. Recently, in f2), the authors con- 
nected mean bit rate in time varying AI/AI/1 queue with moments of the integral Q{u) du in a 
corresponding standard M/M/1 system. 

However, applications of such integrals go beyond queueing systems. Let {S{t),t > 0} be a standard 
risk process. Integrals l{s(^u)<o}S{u) du, where T is deterministic (i.e. integrated negative part 
of the risk process), are suggested in |8l as possible risk measures. Further extensions are given in 
a multivariate setting. Furthermore, as in IflTl , integrals JJ^ Q{u) du in Geo / Geo / 1 queue and (as a 
limit) in M/M/1 system have particular interpretation in compact percolations. Last but not least, if 
X is Levy process, integrals cxp(— du have applications in financial mathematics, see Il25l . 
Another applications are coming from the actuarial science, where very often regulated processes are 
considered and integral functionals from a regulation random mechanism are investigated. 



Consider a stable GI/GI/l queue. Denote by {T,Ti,i > 0} and {S,Si,i > 0} two stationary i.i.d. 
and mutually independent sequences of interarrival and service times, respectively. Let Xt = 1/E[T], 
Xs = l/E[S] and p = Xt/Xs- Let {Q{t),t > 0} be a stationary queueing process and r = mi{t > : 
Q(t) = 0} the corresponding busy period. We shall assume that the distribution F of service time S 
is subexponential (denotes as F £ S). The distribution G is subexponential when 



where G*^ is the convolution of G with itself and G denotes the tail distribution given by G{x) = 



Heuristically, the large area Q[u) du is realized by a customer with large service time 5*0, say, who 
arrives at the very beginning of the busy period and blocks the server. During that time 5*0, ac- 
cording to the Law of Large Numbers (LLN), approximately AtS'o customers arrive and the queue 
length process increases linearly. Hence the area under queueing process before reaching maximum 
is asymptotically equivalent to the area of a triangle: 2l7^(^)^' where Q{t) is maximum of the queue 
length process over the busy period (see Figure 1). After passing the maximum, the queue length 
process behaves according to the LLN, decreasing almost linearly to with the slope Xs — Xt- Thus, 
the area under the queueing process after reaching the maximum is equivalent to 2(As-At) ^^'^-^^ ^^'^ 
Q{t) is equivalent to p{Xx — Xs)t. This heuristic leads to the following conjecture. 

Conjecture 2.1. IfF e S, then 



(1) 




2. Subexponential asymptotics 



lim G*^{x)/G{x) = 2, 



X 



1-G{x). 



(2) 
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Figure 1 . A typical behaviour of a heavy tailed queueing process. 



Similarly, for the workload process the heuristic is as follows. The most likely way for the area to be 
large is that one early big service time occurs and apart from this, everything in the cycle develops 
normally. Using LLN and ignoring random fluctuations, this leads to the conclusion that the work- 
load goes to zero with negative rate — (1 — p). Thus the area exceeds level x iff the area of the triangle 
with the sides t(1 — p) and r is greater than x, hence when 

ir2(l-p) >a;, 

which suggests the following conjecture. 
Conjecture 2.2. IfFeS, then 

(3) p( [ W{u) du>x] p(t> 



2x 

The statements ^ and ^ were proven in l22l and 13/ respectively, under regularly varying assump- 
tion of the service time, that is F(x) = x~°'L{x), where a > 1 and L is slowly varying at infinity. 
Furthermore, in Il22l one needs additionally that 

^ ' t^oo F{t) 

holds with > 0. 

The tail behaviour of the busy period r can be identified in terms of F for a large subclass of S: 

(5) P(r > x) ^ -^F((l - p)x), 

1 - P 

see ||5l, mi and HO). In particular, using (|5) in the regularly varying case, together with Q and (O 
yields exact asymptotics for area under queueing process and workload, respectively (see l22| and 

my 
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3. Light-tailed asymptotics 

As in Section |2] we consider a stable GljGIjX queue. Here, we assume that the service time S is 
light-tailed, that is there exists > such that E[exp(6'S')] < oo. 

Under the above assumptions for the queueing process, we have the following open problem: 
Open Problem 3.1. find exact asymptotics of 

P (^j^ Q{u) du > 

We suppose that 

P Q{u) du>x^ Cx-^^'^ exp(-'(/'Vx) 

for some constants and C. 

We suggest above asymptotics, believing that it should be the same like for M/M/l queue, which 
was found in IT4l under two conjectures on p. 391 and it is in the following form: 

(6) P Q{u) du > ^ l^^x-^l^ exp(-VV^), 

where 

^ = 2v/-2(l-p) + (l + p)logp. 

In the proof authors used the Laplace transform method. We are not aware of any probabilistic proof 
of this result and we do not know if Conjectures 1 and 2 in | [T4| hold true. 

Unfortunately, we have not managed to produce any heuristic for this result either. The idea of the 
piecewise linear most likely trajectory seems to produce wrong expression. In particular, define the 
new probability measure P: 

where J"„ = cr(Ti, 5*1, . . . , r„, 5„) and 7 solves the equation E[cxp(7(T - S))\ = 1. Let p = |||| > 
1. Consider the most likely path coming from large deviation theory for large cycle maxima, that 
is trajectory that develops along the line with the slope p — 1 and after getting maximum behaves 
'normally', that is goes to zero linearly with negative rate —(1 — p). This trajectory produces wrong 
asymptotics for MjMjX queue since by Kyprianou |[T9l we have then: 



P 



[j^^ Qiu) du >x^^p{r> - Cx-3/4e-^\/^v^ 



for some constant C and 7 = (1 — y//5)^/i. The explanation might come from papers | [T2ll26ll27) . where 
the optimal trajectories in a sense of large deviation theory for the mean value of a reflected random 
walk Wn are considered. The relationship between area under the queue length and sum of Wn is 
clear when observing queue length process at arrival and departure epochs. This papers suggest 
that optimal path, though still concave in general, might not be piecewise linear. In fact, this papers 
suggest that the optimal trajectory q{i) of the queueing process for the large value of the area on the 
cycle should solve the following equation: 

(8) yi(^q(^t)]=X*{T-t), 



dt 

where A* > and / is a rate function for the increment process. Still, solving (|8]l explicite seems to 
be a difficult task and transferring it into finding the asymptotic tail distribution of r and the integral 
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even more cumbersome. This problem shows also the need of simulations of the area Q{u) du to 
check if (|8]l indeed produces non-piecewise-linear optimal path. In general, we are not aware of any 
simulations results concerning area under a random process. 
Similar open problem one can pose for P (JJ^ W{u) du> x). 



4. Extensions 

4.1. Transient case. Let us consider the M/M/1 queue with p = 1. In lllTl the author established the 
following asymptotics for the critical case p = 1: 

where ^/jq and il^i are explicit constants. 

Open Problem 4.1. Consider a transient GljGIjX queue such that F{x) is regularly varying. What is the 
tail asymptotics of Q{u) du and JJ" W{u) du ? 

As mentioned above in the Introduction, this type of questions should have cormections with heavy- 
tailed critical percolations. 

4.2. Utility functions, discounting. Consider again a stable G//G//1 queue. Assume that / is a de- 
terministic function. We are interested in the tail behaviour of /J^ f{W(u)) du. In an insurance context, 
/(•) may play a role of the utility function. Of course, one can formulate the corresponding problem 
for the queuing process, where / gives the costs of maintaining the system. 

The heuristic given in the subexponential case in Section |2] suggests that studying the asymptotics of 
latter integral is equivalent to study the behaviour of 

when X becomes large. For example, if one considers f{x) = x^ , /c > 0, it leads to the following 
conjecture. 

Conjecture 4.1. Assume that F <E S and f{x) = x^. Then, 

However, this is not clear for us what should be expected for a very general function /. 

Now, assume that 9 > and consider the process X{t) = exp{—0t)f{W{t)), t > 0. Similarly, the 
corresponding integral can be interpreted in an insurance context: exp(— 0i) is the discounting factor 
and /(•) is the utility function. 
The subexponential heuristic leads us again to 

p(^J^ cxp{-9u)f{{T-u){l- p))du> x^ . 
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For further heuristic, consider f{x) = x*^. Then 

cxp{-0u)f((T - u){l - p)) du = (1 - p)'' I exp(-6lM)(T - u)" du 

= {1 - p)'' exp{-9T) [ exp{9s)s''ds 
Jo 



1 - p)" exp(-6iT) 



■ exp(6'r)T'' 



/ cxp{es)s''-^ ds 
Jo 



Integrating further by parts the last integral we see that the leading term will be 0(cxp(6't)t'^ ^ ). This 
leads to the following conjecture. 

Conjecture 4.2. Assume that F <E S and f{x) ~ x'^. Then, 

p( [ cxp{~du)f{W{u))du>x] ^piT> 



This conjecture should be compared with Conjecture 14.11 showing how discounting leads to the 
change in the asymptotic behaviour. 

Let us consider now a light tailed case. The large area for the integral W{u) du is strictly connected 
to a large W{t). It does not seem to be the case if one considers cxp{—du)f{W{u)) du, i.e. large 
values of W{u) may be killed by the discount factor. 

Open Problem 4.2. Under the conditions of Open Problem \3l] find the asymptotics for 

P exp{-eu)f{W{u)) du> x^ . 

4.3. Finite horizon. All the conjectures and open problems above may be formulated in case of , 
where T is finite and deterministic. In particular, if S{t) = v + ct — J^i'Ji^ ^i' i 5^ 0, is the classical risk 
insurance process, then moments of 

fT 

S{u)l{s{u)<o} du 



are proposed in |8] as particular risk measures. In that paper the authors established the asymp- 
totics for the expected value of the integral, when T is fixed and the initial capital v becomes large. 
Clearly, this problem can be re-formulated for the dual workload process. Furthermore, all the above 
mentioned extensions (discounting, utility functions) may be considered. 



4.4. Multivariate case. As in Section |431 let us consider insurance context. We assume that a com- 
pany has two lines of business and customers arrive according to a renewal process N{t), t > 0, with 
generic interarrival time T. The ith customer has a claim (S*;,!, 5i_2)- For example, a car accident may 
cause a claim for driving and liability insurance (see for example |3l, ||3 and references therein). If 
two lines of business are considered separately, we are interested in the tail behaviour of 



Si{u)l[Si{u)<Q}du, J S2{u)l[S2(u)<a}du j , 

where 5*1 and 5*2 are risk process associated with the corresponding lines of business. One can con- 
sider the tail behaviour when T is fixed and a vector of initial capital {xi,x2) becomes large (in a 
particular sense, usually dependence is linear, that is xi = ax2 for some fixed a > 0). 
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In queuing theory one can consider e.g. parallel queues, where customers are coming into the system 
according to renewal process N{t) and each service time (we assume here that all of them are i.i.d.) is 
proportionally divided into two servers (see e.g. |23| and Il24l ). In this case, Si,i = 65^2 for some fixed 
b. Let Wi{t) and Qi{t) will be the workload and queueing process on ith server (i = 1,2). The busy 
period t is understood here as the minimum of the busy periods on both servers. We could analyze 
the tail behaviour of the following "tail" probability: 

P Wi{u) du > X and J W2{u) du > axj as a; oo, 

under assumption that distribution F of Si.i is subexponential or light-tailed. Similar considerations 
could be analyzed for the bivariate occupation process. 



4.5. Regulated processes. Apart of original process X{t) in ((T]|, one can consider regulated process: 

U{t) = X{t)-L{t), 

where L{t) is left-continuous, increasing process, adapted with respect to the natural filtration of X. 
For example, if X is a spectrally negative Levy risk process, then L{t) might be cumulative dividends 
payment up to time t payed according to some strategy. The most often used strategy is so-called 
barrier strategy in which all surpluses above a given level a are transferred (possibly subject to a 
discount rate) to a beneficiary. In this case L{t) = a V sup^<j ^i-'') ^ « ^rid for given utility function 
we are interested in in the following random variable: 

(10) I fir) = r e-'y{L{u))du, 

Jo 

(possibly with 6 = 0), where t = inf{t > : U{t) < 0} is a ruin time. 

We can also analyze discounted cumulative dividends payed up to ruin time: 



(11) / e-''"dL(u), 

Jo 

which is equivalent to finding the tail asymptotics of integral: 







where {Zi{t), Z2{t)) is a bivariate subordinator and AZ2{u) = Z2{u) — Z2{u—) is a size of jump of Z2 at 
time u (see |TH4l|20H29l for details and other references). In the case if a = 00 we end up with integral 
from exponential function of Levy process, see e.g. (51[71[10j|25| for this kind of functional. Note that 
in particular / = g^sziit) solves the following equation: / = /J" e~^^i(*) du + e~^^^^'^^I and 
analyzing / is strongly related then with properties of e"^^^*^"^) (see e.g. IITSllTSl ). 
With this problem there is also related integral with respect to the general Levy process: 

cT 

Y{u)dX{u) 







for appropriate predictable integrand Y and fixed T (see IITSl for multivariate regularly varying set- 
ting). 

In queueing systems with single server and finite capacity a, L{t) corresponds then to cumulative 
lost of work. If only proportion of information is lost, then we can consider e.g. L{t) = cl{t/(t)>a} 
for some c > 0. In this case formally so-called refracted process U solves the following stochastic 
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equation: 

dU{t) = dX{t) - cl{a(t)>a} dt; 
see IIZTl . In this case we can also try to find asymptotic tail of flOt and i fTTI l. 
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